Symmetry, Integrability and Geometry: Methods and Applications 



SIGIVIA 5 (2009), 069, 17 pages 



Non- Hermit ian Quantum Systems 

and Time-Optimal Quantum Evolution* 

Alexander I. NESTEROV 

Departamento de Fisica, CUCEI, Universidad de Guadalajara, 
Av. Revolucion 1500, Guadalajara, GP 444^0, Jalisco, Mexico 
E-mail: nesterov@cencar.udg.mx 

Received November 17, 2008, in final form June 23, 2009; Published online July 07, 2009 
doi:10.3842/SIGMA.2009.069 

Abstract. Recently, Bender et al. have considered the quantum brachistochrone problem 
for the non-Hermitian PT-symmetric quantum system and have shown that the optimal time 
evolution required to transform a given initial state into a specific final state j'i/'/) can be 
made arbitrarily small. Additionally, it has been shown that finding the shortest possible 
time requires only the solution of the two-dimensional problem for the quantum system 
governed by the effective Hamiltonian acting in the subspace spanned by j'i/'i) and | ■(/'/)• In 
this paper, we study a similar problem for the generic non-Hermitian Hamiltonian, focusing 
our attention on the geometric aspects of the problem. 

Key words: non-Hermitian quantum systems; quantum brachistochrone problem 
2000 Mathematics Subject Classification: 81S10; 81V99; 53Z05 



1 Introduction 

In view of recent results on optimal quantum evolution and its possible relation to quantum 
computation and quantum information processing, there has been increasing interest in the 
quantum brachistochrone problem. The problem consists of finding the optimal time evolution r 
to evolve a given initial state into a certain final state {ipf) under a given set of constraints [1, 
2, 3, 4, 5]. It is known that for the Hermitian Hamiltonian, r has a nonzero lower bound. 
However, Bender et al. [2] have recently shown that for non-Hermitian "PT-symmetric quantum 
systems, the answer is quite different and that the time evolution r can be made arbitrarily 
small, despite the fact that the eigenvalue constraint is held fixed and identical to that for the 
corresponding Hermitian system. The mechanism described in [2] resembles the wormhole effect 
in general relativity and has generated discussion in the literature [6, 7, 8, 9, 10, 11, 12, 13, 14]. 

Non-Hermitian Hamiltonians emerge in physics in different ways. For instance, the applica- 
tion for description of dissipative systems is well known from the classical works by Weisskopf 
and Wigner on metastable states [15, 16, 17]. It was demonstrated that the evolution of the 
quantum system, being initially in the metastable state il^{0), can be described by the effective 
non-Hermitian Hamiltonian ffeff as follows: V'(O) ~^ V'(^) = e~*^°"*^(0)-|- decay products. Re- 
cently, it has been shown how the non-Hermitian Hamiltonian appears in the framework of the 
quantum jump approach to open systems [18, 19]. Other examples include complex refractive 
indices in optics, complex potentials describing the scattering of electrons, atom diffraction by 
light, line widths of unstable lasers, etc. For the following, it is essential that non-Hermitian 
physics differs drastically from the conventional physics in the presence of the so-called excep- 

*This paper is a contribution to the Proceedings of the Vllth Workshop "Quantum Physics with Non- 
Hermitian Operators" (June 29 - July 11, 2008, Benasque, Spain). The full collection is available at 
http://www.emis.de/journals/SIGMA/PHHQP2008.html 
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tional points, where the eigenvalues and eigenvectors coalesce, even if the non-Hermiticity is 
regarded as a perturbation [20] . 

In the Hermitian quantum mechanics, the optimal time evolution problem implies finding 
the transformation j-f/jj) — > e~*^*|V'j) that provides the shortest time t = t under a given set of 
constraints [1]. The generic non-Hermitian quantum brachistochrone problem poses the same 
question, with the exception that now the evolution of the system is described by the non- 
Hermitian Hamiltonian, which is not necessarily PT-symmetric [2, 6, 12]. 

Recently, this problem has been studied by Assis and Fring for a dissipative quantum system 
governed by a symmetric non-Hermitian Hamiltonian [6]. It has been shown that the passage 
time required to transform a given initial state to the orthogonal final state can be arbitrarily 
small. The obtained effect, being similar to the one discovered by Bender et al. [2], is related 
to the non-Hermitian nature of the system rather than to its PT-symmetry. 

In this paper, we address the non-Hermitian quantum brachistochrone problem, considering 
the generic non-Hermitian Hamiltonian, and focus our attention on the geometric aspects of the 
problem. In Section 2 we first briefly summarize the main properties of non-Hermitian quantum 
systems. We discuss the non-Hermitian quantum brachistochrone problem and the associated 
two-dimensional effective Hamiltonian, acting in the two-dimensional space spanned by the 
initial and final states. In Section 3 we discuss the quantum brachistochrone in the vicinity of 
the exceptional point. In Section 4 we consider the non-Hermitian quantum brachistochrone 
problem in the context of the Fubini-Study metric on the complex Bloch sphere. In Section 5 
the results and open problems are discussed. 



2 Non-Hermitian quantum brachistochrone problem 

In the quantum brachistochrone problem, finding the shortest possible time requires only the 
solution of the two-dimensional problem, namely, finding the optimal time evolution for the 
quantum system governed by the effective Hamiltonian acting in the subspace spanned by \ipi) 
and IV-/) [2, 3, 4]. 

Before proceeding further, we outline some background information on non-Hermitian quan- 
tum systems. Let an adjoint pair {\tp(t)), \tp{t))} be a solution to the Schrodinger equation and 
its adjoint equation 

which can be recast in the form 

For Afc being the eigenvalues of H, we denote by |'i/'/c) and {ipkl the corresponding right and left 
eigenvectors: H\tpk) = AfclV'fc)) {ipk\H = Xkitpkl- The systems of both left and right eigenvectors 
form bi-orthogonal basis [21] 

Let the set \ipo)^ (V'ol} form the bi-orthonormal basis of the two-dimensional sub- 

space spanned by the initial state l^/^j) and the final state \ipf)'- 
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Using this basis, we can write the final state {ipj) as 



|^;)=cos||^,)+e^^sin||V'o), (2.2) 



where a, (3 are complex angles, and similarly. 



(V^/l = cos + e *'^sin|(V;o|. 

Let \ip{t)) = ui{t)\'ipi) + U2{t)\tpQ) and {ip{t)\ = ui{ilji\+U2{Tpo\ be solutions to the Schrodinger 
equation and its adjoint equation (2.1), respectively. Then one can see that the vectors \u(t)) = 

^^1^1 ^ {u{t)\ = {ui{t) , U2{t)) satisfy 

ig-^Ht)) = H,if\u{t)), = {Ht)\HeH, (2.3) 

where the effective two-dimensional Hamiltonian f^eff reads: 
1 f Xo + Z X-iY 



^"^^^ 2yX + iY Xo-Z 
Here, we set 

Ao = {HHm + {MH\ipo), X = {A\H\iPo) + {MHm, 
Y = i{{A\H\iJo) - (HHrn), Z = {^P,\H\A) - {i'o\H\i^o). 

Further, it is convenient to express the effective Hamiltonian ifgff in terms of the Pauli 
matrices: 

H,if = ^l + ^n-a, (2.4) 

where 1 denotes the identity operator and = (X, Y, Z) . 

It is then easy to show that the complex Bloch vector defined as n(t) = {u{t)\a\u{t)) , where cr 
are the Pauli matrices, satisfies the complex Bloch equation 

which is equivalent to the Schrodinger equation (2.3) (see, e.g., [22, 23]). In the explicit form, 
the components of the Bloch vector are given by 

ni = uiU2 + U2U1, n2 = i{uiU2 - U2U1) , = uiui - U2U2. (2.5) 

The vector n(i) = (ni(t), n2(t), n3(t)), being a complex unit vector, traces out a trajectory 
on the complex 2-dimensional sphere , and the latter can be considered as the quantum phase 
space for the non-Hermitian two-level quantum system. From equation (2.5), we find that 
rii = (0, 0, 1) corresponds to the initial state 

Analysis of the eigenvalue problem He^f\u±) = A|n)-|-, (n±|i?eff = -^±(^^±1) yields X± = 
(Ao ± i?)/2, where R = ^JX'^ + Y'^ + Z"^. The right and left eigenvectors are found to be 

l^+^ = (e^'Cin| )' -2 

-e-'^sinf \ 

cos^ r ' V " "2' --'2 
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where 



R+Z . 9 R-Z 



2 y 2R 2 y 2R 

e- = ^±^, e- = ^^, (2.7) 



and 

X = ii sin cos 99, Y = RsinOsimp, Z = Rcos9, (2-8) 

(0, 93) being the complex spherical coordinates. 

The coalescence of eigenvalues A+ and A_ occurs when + Y'^ + Z"^ = 0. There are two 
cases. The first one, defined hy 9 = 0, ip = 0, yields two linearly independent eigenvectors. The 
related degeneracy is known as the diabolic point, and we obtain 

ly (n+| = (l,0), |n_) = (°), (n_| = (0,l). 

The second case is characterized by the coalescence of eigenvalues and the merging of the eigen- 
vectors. The degeneracy point is known as the exceptional point, and we obtain: = e^'^\u-) 
and (u+l = e~^'^{u-\, where k G C is a complex phase. 

Let us assume that the exceptional point is given by Rq = {Xq,Yo, Zq). Then, if Zq / 0, 
using equations (2.6)-(2.8), we obtain 

tan^ = ±z, e^^^" ^° + 



2 ' Xo-iYo' 

thus, at the exceptional point 90 ±00. If Zq = 0, we obtain Xq = HYq. This implies that 
^0 = 7r/2, and Qip — > ±00 at the exceptional point. 

Taking the Hamiltonian of equation (2.4), we find the solution of the Schrodinger equa- 
tion (2.3), satisfying |V'(0)) = IV'i)) as 

\m) = Ci{t)e-'^ot/2^i;i) + C2(t)e-^^«*/Vo), (2.9) 

where 

Ci(t) = cosy -icos6'siny, C2(t) = -ie*^ sin 0sin ^ , (2.10) 



and we denote n = R = VX'^ + Y'^ + Z"^. 

The solution of the adjoint Schrodinger equation with the wave function {ip{t)\ written as 

m)\ = Ci(t)e^^»*/2^^^| + C2{t)e''^'^'/^{H (2.11) 
is given by 

C'i(t) = cos Y -hicos6i siny, ^2(0 = «e"'^ sin6l sin y . (2.12) 



Applying equation (2.2), we can write \ip{t)) as 

-i/3 

a 
2 



m)) = (Ci{t) - e-'^cot "C2{t)) e-*^°*/V.) + ^C'2(t)e-*^«*/V/>- (2.13) 
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Hence, the initial state {ipi) evolves into the final state {ipf) in the time t = t when 

Ci{t) - C2{T)e-'^ cot 1 = 0, C2(r) = e^^sin |. (2.14) 

To find the time evolution r, we first solve the equations (2.14) for e**^'^"'^). The computation 
yields 

cos^ 6 — s\v? ^ 



sm sm ^ 



— cos 6 cos § ^ A /cos^ — sin^ ^ 

e"^(^-'^) = ^. (2.16) 

sm sm 2 

It should be noted that the solution related to the lower sign can be obtained from the solu- 
tion corresponding to the upper sign by changing the parameter a as follows: a — > a -|- 2ti. 
This implies the change of the total sign in the final state of the quantum-mechanical system: 
\iljf) —\ipf)- Thus, without loss of generality, one can consider only one sign (upper or lower) 
in equations (2.15), (2.16). Further, for definiteness, we will choose the upper sign. Then, 
substituting the result into (2.14), we get 

r2r i sin^ f 
tan — = , ^ . (2.17) 

cos I -V / cos^ 9 — sin^ ^ — cos 9 
From here, the time evolution r is found to be 



2 / isin^f 



arctan 



(2.18) 



^ y cos ^ ^ cos^ 9 — sin^ ^ — cos 

In addition, since r should be a real positive function, the following restriction must be imposed: 

(i sin^ — 
^ I . (2-19) 
cos tta/cos^ 9 — sin^ § — cos ( 



Now, the generic problem is to select a final vector \ipf), choosing the parameters a and /3. 
Next, we must find the conditions that should be imposed on the parameters {9,ip) to yield 
the smallest time r required to evolve the state into the state \^f) under a given set of 
constraints with the Hamiltonian 

^•^^ = 2 V Ao j + 2 V e^^sin^ -cos9 ) ' ^^'^^^ 

In what follows, we assume the eigenvalue constraint to be imposed as |A+ — A-| = \^\^. This 
implies that the argument of il is determined by equation (2.19). 

Further study of the critical points shows that there is no solution with a finite value of |0| 
yielding the minimum of the time evolution [24]. Moreover, as follows from equation (2.18), in 
the limit |^| oo (|9^| oo), the time evolution behaves as r ~ |2/r2cos0| 0. Thus, for 
a quantum-mechanical system governed by a non-Hermitian Hamiltonian, the time evolution r 



^In contrast to the Hermitian quantum systems, here, there are more choices for imposing the eigenvalue 
constraint. For instance, instead of | A+ — A- 1 = |r2| being held fixed, one can require A+ — A- = = const [24]. 
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can be taken to be arbitrarily small. In addition, since for any finite value of |^|, the minimum 
of r does not exist, the non-Hermitian Hamiltonian cannot be optimized. 

However, for a quantum-mechanical system governed by a Hermitian Hamiltonian, the latter 
can be optimized. Indeed, in the Hermitian case, Q(p = = 0; hence, the saddle point 6 = tt/2 
becomes the point of a local minimum (see Fig. 1). It then follows from (2.20) that 

_i(Xo \ n( e-^^ 

~ 2 V Ao y 2 V e*^ 
and from equations (2.13)-(2.16), we obtain 

\m) = e"*^"*/' (^cos — - cot - sm — J jV'i) + — ^ — jV'/). 

This fully agrees with the results obtained by Carlini et al. and Brody and Hook [1, 4]. 

We now consider some illustrative examples. Of a special interest is the case when the 
complex Bloch vector n, entering in the non-Hermitian Hamiltonian (2.4), is orthogonal to the 
plane spanned by rij and n f . Without loss of generality, one can choose n = rii x Uf / sin a and 
set rij = (0,0, 1). This yields 9 = Tr/2 and from equation (2.15), we obtain e^^ = ie^^ . Applying 
these relations, we get 

Then, using equation (2.17), we obtain the optimal time evolution as r = |a/r2|, and from 
equation (2.19), we get Q = |r2|e*'^''^°. This example has a simple geometric interpretation. 
As can be easily shown, the final state nj \s obtained from the initial state rij by rotating the 
complex Bloch sphere S"^ through the complex angle a around the axes defined by the vector n. 
It should be noted that for SJa = 9/3 = 0, we have Q'il = 0, and the Hamiltonian (2.21) coincides 
with the optimal Hamiltonian, yielding the shortest time evolution for the unitary evolution. 
As we will show in the following sections, Tm gives the upper bound of the time evolution for 
the generic non-Hermitian Hamiltonian. 

The other interesting example is when the initial and final states are orthogonal to each other 
{a = vr). The smallest time Tp required to evolve from a given initial state iV'i) to the orthogonal 
final state is called the passage time [3, 4]. Inserting a = vr into equation 2.18, we obtain 

^ In 



9. z + n 

Let us consider some limiting cases, starting with the saddle point Z = (see Fig. 1). 
Computation yields Tp = tt/\0,\. This can be interpreted as the Hermitian limit of the non- 
Hermitian system. As can be seen in Fig. 1, the evolution of the non-Hermitian quantum 
system is faster than that of the corresponding Hermitian system, satisfying the same eigenvalue 
constraint. Moreover, the passage time Tp — > while \Z\ — > oo. 

Next, we find that — > oo at the points Z = ±0 (see Fig. 1). To understand this result, we 
note that Z = 0, implies 9 = and that Z = 0, yields 9 = tt. For both cases, 9 = and 9 = n, 
the wave function (2.13) becomes 

Thus, only phase and amplitude of the wave function are changed; otherwise, the system in the 
same initial state remains intact. This explains the divergence of the passage time at the points 
Z = zizQ. This is easy to understand by looking at the example of the half-spin particle in the 
uniform magnetic field, when the spin is eqnarrayed with the magnetic field. 
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Figure 1. Plot of the passage time Tp versus diZ and Left panel: diTl = 1, Sfi = 0.1. Right panel: 
SRfi — 1, — 10. As can be seen from the plots, the singularity occurs at the points Z = ±il and 
Tp = 7r/|il| at the saddle point Z — 0. 



3 Quantum evolution in the vicinity of the exceptional point 

In the presence of exceptional points characterized by coalescence of eigenvalues and correspon- 
ding eigenvectors, non-Hermitian physics differs essentially from Hermitian physics [20, 25, 26]. 
Since for the Hermitian operators, the coalescence of eigenvalues results in different eigenvectors, 
exceptional points do not exist for Hermitian systems. The related degeneracies, referred to as 
'conical intersections', are known also as 'diabolic points'. In the context of Berry phase, the 
diabolic point is associated with the 'fictitious magnetic monopole' located at the diabolic point 
[27, 28]. In turn, the exceptional points are associated with the 'fictitious complex magnetic 
monopoles' [29]. 

For a two-dimensional system (2.4), the eigenvalues coalesce when $7 = 0. Writing the 
complex vector ft as ft = r — iS, where r and S are real, one can recast the effective Hamiltonian 
of equation (2.4) in the form 

H,f = ^l + \r-a-'-6-cT. (3.1) 

As can be seen, the degeneracy points are defined by the following equation: 

-(^2 -2ir(5cos7 = 0, (3.2) 

where 7 denotes the angle between 6 and r = (x,y,z). Furthermore, it is convenient to use 
cylindrical coordinates {p,z,(p), in which 0, = p"^ + z"^ — 5"^ — 2iz6, and equation (3.2) reads: 
p'^ + z'^ — 6'^ — 2iz6 = 0. There are two solutions to this equation. The first one is p = z = 6 = 0, 
and the corresponding degeneracy point located at the origin of coordinates is the diabolic 
point. The second solution, given by z = and (? — b'^ = 0, defines the exceptional point. 
Generally speaking, the exceptional "point" , being realized as a cone in the four-parameter space 
{x,y,z,5), is not a point in the convenient sense, and all possible one(zero)-dimensional cases 
can be obtained by the conic sections. For instance, for a fixed 6, we obtain the exceptional point 
as a circle of radius \5\, lying in the plane z = 0. It should be noted that in the recent literature, 
the term "exceptional point" is applied not only to the particular case when the exceptional 
"point" is indeed a point but to the general case when the non-Hermitian degeneracy is realized 
as a submanifold in the parameter space [30, 31, 32, 33]. In the following, we will stick to this 
more general interpretation of the exceptional point. 
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As follows from equations (2.9)-(2.13), at the exceptional point, the quantum evolution is 
described by 

Without loss of generality, we further assume 6 > 0. The computation of the time evolution at 
the exceptional point then yields 



2 

^=5 



a 

1-cos- 



In what follows, imposing the eigenvalue constraint as \Q\ = const, we restrict ourselves 
by consideration of the orthogonal initial and final states Substituting a = tt into 

equation (2.13) and taking into account equations (2.15), (2.16), we obtain 

|V'(t))= (^cosy -icos^siny^ e-^^o^lV'i) +sin0sinye-^^o*|'0/). (3.3) 

Now recalling that cos^ = (z — i6)/0, and setting z = 0, we observe that in the vicinity of the 
exceptional point, there are two different regimes, depending on the relation between p and S. 

For p > 6, we have 0, = (?■ — > 0, and inserting cos0 = (z — iS)IVl into equation (3.3), 
we obtain 

m)) = (cos ^ - ^ sin ^) e-^o.|^^) + A ^e^^o*,^^), 

where JIq = — (^^|^^^ denotes the Rabi frequency. If, in addition, Q'Ao = 0, then the eigenvalues 
of the non-Hermitian Hamiltonian (3.1) are real, and we obtain the "PT-symmetric Hamiltonian 
widely discussed in the recent literature in connection with the quantum brachistochrone problem 
[2, 5, 7, 8, 9, 10, 11, 12, 34, 35]. Then, applying (2.18), we get (see Fig. 2) 



2 



Tp = ^ arctan ( ^ ) . (3.4) 



As can be shown, the passage time Tp is bounded above and below as follows: 2j \J + < 
Tp < min{7r/r2o) 2/(5}. This improves the estimation of the passage time obtained in [2]. 

As follows from equation (3.4), for ^ 6, the passage time has the same behavior as for 
the equivalent Hermitian system: Tp ~ it/0,q, and if 0,q <^ 5, one obtains Tp ^ 2/6. At the 
exceptional point = (and p = 6): 



From here, we get Tp = 2/5 and ^ if (5 — > oo. Thus, in contrast to the Hermitian systems, 
the passage time may be arbitrarily small; this, of course, implies that p, 6 becomes large. This 
is possible due to the hyperbolic nature of considered as function of p, and 5 and agrees with 
the results obtained in Bender et al. [2]. As a result, the equivalent Hermitian system in the 
limit Oo — > yields Tp = tt/Qq oo. 

Now assuming p < 6, we have O = iO,o and 
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Figure 2. Plot of passage time Tp vs. x — S {fto — 1). Left panel: — (2/rio) arctan(rio/^) (solid 
green line), t = + (dotted blue line) and t = 2/6 (dashed red line). Right panel: Tp = 

(2/rio) tanh^^(rio/^) (solid red line) and r = 2/S (dashed green line). It is seen that in both cases, Tp is 
asymptotically approximated by t = 2/d. 



Computation of the passage time yields 

- - 2 tanh-if^V (3.5) 



^^0 \ ^ , 

As seen in Fig. 2, the passage time is bounded below as follows: Tp > 2/5, and it follows from 
equation (3.5) that ~ 2/5 in the limit 0,q ^ 5. At the exceptional point, we obtain the same 
result as above: Tp = 2/6 and 

l^(t))=(l-|) e-^^°*|V'.) + |e-^^"*|V'/). 

As an illustrative example, we consider a two-level dissipative system driven by a periodic 
electromagnetic field E(i) = 3f?(£(t) exp(iz/t)). In the rotating wave approximation, after remo- 
ving the explicit time dependence of the Hamiltonian and the average effect of the decay terms, 
the Schrodinger equation reads [36, 37] 

ui\ ^If -iX + A-i6 2V* (o fi^ 

U2 J 2 V 2y -iX-A + i6 J {u2 J ' ^ ^ 

where A = {ja + 7b)/2, with 7^, 7^ being decay rates for upper and lower levels, respectively, 
A = Ea — Eh — v, 6 = (7a — 7;,) /2, V = ■ S , and /x is the dipole matrix element. 

The choice S{t) = £oexp{iu;t) yields V{t) = Vq exp{iujt) , where Vo = fj, ■ So, and we assume 
further that Vq > 0. The solution of equation (3.6) with this choice of S can be written as 

\u{t)) = Ci(t)e-*('"-^^)*/>t) + C2(t)e^(^+^^)*/>^), 
where \u^) ~ ^ q ^ ^^"^ 1^1^ ~ ( 1 ) ' up/down states, respectively. 



\C{t)) 



Ci(t) \ _ / cosf -icos^sinf -isin^sinf \ / Ci{0) 
C2{t) ) ~ \ -i sine sin f cos ^ + icos^sin f J V Csio) 



and we set cos 6* = (A - a; - i6)/n, n = {p"^ + {A - lv - iSf)^/'^, p = 2Vb. 
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As can be shown, the state \C{t)) satisfies the Schrodinger equation 

AC) 



dt 



Hr\C), 



written in a co-rotating reference frame, where the Hamiltonian Hr takes the form 



Hr 



n 

2 



cos ( 
sin ( 



sm U 
- cos ( 



Let \u{t)) be the solution of equation (3.6) with the initial state at t = taken as \u^) and 
the final state of the system at a later time t being \u-i) or Then following Baker [38], we 
compute the probability (-Pit) system is in the state \ui) (|nj)), respectively: 



TT 



cos 



I COS t> sm 



-At 



Pi 



it 



sm t) sm 



-xt 



For the resonance frequencies, w = A, we have 0, = {p^ — 5^)^/^, and cosO = —i5/^Q. There 
are two different regimes, depending on the relation between p and 8. For p > 5, we have 
a coherent tunneling process 



P^ 



TT 



cos 



6 not 

— sm 

no 2 



Pi 



not 



TT 



-xt P ■ 2 "0 

e ^sm — 

nl 2 



(3.7) 



where 
incoherent tunneling 



p2 _ j2|i/2 jg22otes the Rabi frequency. On the other hand, for p < 6, there is 



TT 



-xt 



cosh 



not 



6 , not 
n-o''""^^ 



Pi 



not 



TT 



{3i 



At the exceptional point, JIq = 0, and both regimes yield 



TT 



6t 



-xt 



Pi 



TT 



6t 



-xt 



This is in accordance with the results obtained by Stafford and Barrett and Dietz et al. [39, 
40]. Moreover, the decay behavior predicted by equations (3.7)-(3.8) has been observed in the 
experiment with a dissipative microwave billiard [40]. 

We find that for p > 5, the passage time required to transform the state \ui) into the state |uj) 
is given by Tp = {2/uJo) arctan(u;o/'5)- Similarly, for p < 5, we obtain Tp = {2/ujo) tanh~"^(a;o/<5). 
At the exceptional point, for both regimes, the passage time is found to he Tp = 2/6 (see Fig. 2). 



4 Fubini— Study metric on the complex Bloch sphere 
and the brachistochrone problem 

In this section, we will consider the quantum non-Hermitian brachistochrone problem in the 
context of the geometric approach developed by Anandan and Aharonov [41]. Let \ip{t)) and 
{ip{t)\ satisfy the Schrodinger equation and its adjoint equation, respectively: 

d\m) = H\m), -d{m\ = {m\H, 



Non-Hermitian Quantum Systems and Time-Optimal Quantum Evolution 



11 



where H is a non-Hermitian Hamiltonian, and the standard normalization is held, {ilj{t)\'ilj{t))=l. 
We define the complex energy variance as 



Now, applying the Taylor expansion to {ip^t + dt)) and using equation (4) and its time 
derivative, we obtain 

{'4){t)\'^{t + dt)f = 1 - AE^dt^ + 0{dt^). 

Then, introducing the complex metric as ds^ = 4(1 — {ip{t)\ip{t + dt))"^), we obtain 

ds^ = AAE'^dt^ = Adsls, (4-1) 

where 

d4s = W|(l-P)|#) (4.2) 

is a natural generalization of the Fubini-Study line element to the non-Hermitian quantum 
mechanics, P = \'tp){^l^\ being the projection operator to the state It is easy to show that the 
complex-valued metric (4.2) is gauge invariant, i.e., it does not depend on the particular choice 
of the complex phase defined by the map: {ip) — > e*"!^/'), {tp\ — > e~*"('0|, a € C. We define the 
distance between two given states \tpo) and 1-01) as 

s = '^j \AE{t)\dt, (4.3) 

where the integration is performed along a given curve C connecting \iPq) and 

In the two-dimensional case, the complex-valued Fubini-Study element has a nice geometric 
interpretation. We define a complex distance between two nearby Bloch vectors n(x) and 
n{x + dx) by 

A^(x, X + dx) = 1 — n(x) • n(x -|- dx). 

Then, applying Taylor expansion, 

n(x -I- dx) = n{x) + din{x) dx^ + ^didjn{x) dx'' dx^ + • • • , 
and using n • n = 1, we obtain, up to second-order terms, 

n{x + dx) = 1 — dn{x) ■ dn{x). 
This yields 

ds'^ = dn ■ dn = gij dx' dx^ , (4.4) 
where gij = din ■ dju, and the length of any curve C on S'^ is given by 

L{C) = j Wdn-dn\. 

Denoting n = (sin C cos i^, sin sin z^, cos C), where C,,v we recast (4.4) as 

ds^ = dC^ + sm^ Cdv"^ . (4.5) 
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Finally, using the definition of the complex Bloch vector n = we find that the metric 

on the complex Bloch sphere S'^ can be written as 

ds^ = dn-dn = Adsl^ = 4(dV'|(l - P)\dtij), 

where ds^g = {d'4)\{l — P)\dip) is the Fubini-Study line element. 

Strictly speaking, gij being a complex-valued tensor does not define a proper metric on the 
complex Bloch sphere. However, the advantage of this definition is that, contrary to the Kahler 
metric, the complex "metric" (4.4) is invariant under the gauge transformation e^°'\'ip), 
{ip\ — > e~*"(^/'|, where a G C 

For the Hamiltonian (2.20), the straightforward computation yields 

2AE = nsme, (4.6) 
and, using equations (2.18) and (4.3), we obtain 



L = 2\AE\t 



I ^ sin^ f 

2sin6'arctan 



cos ^ W cos^ 9 — sin^ | — cos ( 



(4.7) 



To compare our results with the PT-symmetric model widely discussed in the literature (see, 
e.g., [6, 7, 8, 9, 10, 11, 12, 13, 14]), we choose 9 = tt/2 + ir] and assume QXq = Qip = to write 
the effective Hamiltonian of equation (2.20) as 

^l/^^o 0\1^ / -is'mhr] e-*^cosh7?\ , . 

''^^ 2 \ Xo J ^ 2 \ e'^coshr] isinhr? J ' ^ > 

Next, we set Aq = rcos7, $7 sinh = rsin7, and ilcoshry = uj. Denoting 6 = rsin7, we 
obtain the Hamiltonian of the system in its conventional form (see, e.g., [2, 5]): 

_lf re-'^ cu e"*^ 
2 \ uj e'"^ re'^ 

where VI = \/ up' — svt? 7 = \JiJ^ — S^. In this particular case, the complex Bloch sphere 
becomes the one-sheeted two-dimensional hyperboloid with the indefinite metric given by [24] 

ds^ = cosh^ p du'^ - dp^, (4.9) 

where —00 < p < 00 and < < 2tt are the inner parameters on H^. It should be noted that 
the interval (4.9) can be obtained from equation (4.5) by the substitution ( 7r/2 + ip, and we 
assume = 0. 

The amount of time r required to evolve the initial state {ipi) into the final state \ipf) can be 
found from equation (2.18) by substituting 9 = -k 12 -|- ir\. The computation yields 

2 / sin^f 
r = - arctan | , (4.10) 

\ cos I -V / sinh^ r] -|- sin^ -|- sinh r/ ^ 
and from equation (4.7), we obtain 

L = 2 cosh r/ arctan 



sin^ I 



cos S^i/sinh^ r\ -|- sin^ % -|- sinh r] 
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To study the spin-flip, we choose the initial state as = |tt|). Then, substituting a = n 
into equation (4.10), we find the time interval 

2 1 2 n 27r 2 17 

Ti = — arctan — ^ — = — arctan — , '^-\ = — arctan — , 
U sinhr/ U o U U 

necessary for the first spin flip from up to down and back. For all values 0, € [0,oo), we 
have 2/6 < < vr/Jl. Thus, the passage time Tp = Tj lies below the Anandan-Aharonov lower 
bound, Tj < 7r/J7 for a spin-flip evolution in a Hermitian system [42], and Tp reaches its minimum 
value Tmin = 2/6 at the exceptional point. In addition, the total time for a spin-flip evolution, 
IT) ~^ I i) ^ I T)i remains invariant: t = Tj -|- T| = 2tt/^1. This is in accordance with the 
results of previous studies [34, 10]. 

Introducing the variable k = aictaii{6/Q) = arctan(sinhr/), we reproduce our results in the 
more familiar form, widely known in the literature (see, e.g., [2, 3, 4, 10]): 

vr — 2k it + 2k 

We see that in the Hermitian limit, r/ — > (5 ^ fi) that implies k ^ 0, the passage time is 
given by Tp = tt/O. In the other limiting case ry — > oo (17 ^ 5), the angle k approaches 7r/2, 
and Tp tends to zero. In terms of variable k, the relation 0, cosh rj = u becomes 17 = a; cos k. 
Furthermore, if the energy constraint — E- = 17 is held fixed, in order to have the passage 
to the limit k — > '7r/2, one must require ^ 17. It then follows from the relation 17^ = u;^ — (5^ = 
const that one must require 5 ^> 17. 

Similar consideration of the distance between the initial and the final states yields 

(1 \ vr — 2k 
— 1 — ) = • (4-11) 
sinh rj J cos k 

In the limit n vr/2, we get Lp — > 2, and in the Hermitian limit, k ^ 0, we have Lp vr. It 
then follows from equation (4.11), that the distance between \u^) and juj), being measured on 
the one- sheeted two-dimensional hyperboloid with the indefinite metric of equation (4.9), is 
bounded as follows: 2 < Lp < tt. 

Inserting 6 = Tr/2 + irj into (4.6), we obtain 2AE = 17 cos??. Then substituting this result 
into equation (4.1), we find that for the non-Hermitian Hamiltonian (4.8), the evolution speed 
V = ds/dt is given hy v = uj = 17 cosh r/. Hence, v > Vg, where = 17 is the speed of the 
geodesic evolution [24]. Similar consideration of the quantum-mechanical system governed by 
the Hermitian Hamiltonian yields v = 17 sin ^, and, obviously, v < Vg. This proves that non- 
Hermitian quantum mechanics can be faster than Hermitian quantum mechanics. Moreover, 
since for any complex angle 9, there exists Oq determined by the equation cos ?R.6q = sinh Q9q 
such that V = |17sin^| > |17| = Vg \Q9\ > \QOq\, this conclusion is applied to an arbitrary 
non-Hermitian Hamiltonian (see Fig. 3). 

Our results are in agreement with those obtained previously by Bender et al. [2]. Howe- 
ver, interpreting the obtained results for r requires care. Indeed, as was pointed out by 
Mostafazadeh [7], to compare the time evolution for the non-Hermitian and Hermitian Hamil- 
tonians and to conclude in which case the evolution is faster, one should not only impose the 
same set of constraints for both cases but also fix the geodesic distance between initial and final 
states. 

In what follows, we consider in detail a spin-flip for the Hamiltonian (4.8) with two types 
of constraints: a) the energy variance AE^ = {ip\H'^\t(;) — {{Tp\H\Tp))'^ held fixed [1], and b) the 
difference between the largest and smallest energy eigenvalues has a fixed value, E^ — E- = 
const [2]. In terms of the parameters (17, a;), the energy constraints are written as follows: 
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Figure 3. Plot of the evolution speed Iwl/I^ij versus x ^ and y = ^9. As can be observed, for any 
e : <m <Tr, there exists the angle Oq such that \v\/\n\ > 1, if |36'| > |96'o|. 



2AE = uj and Ej^ — E- = 0,. Then, using the relation coshr/ = uj/0,, we find that the passage 
time can be recast as 

Tn = 7z arctan =. (4.12) 

As can be observed, under the constraint lu = const, the time evolution is bounded by 2/uj < 
Tp ^ tt/lo. Here, the passage time reaches its minimum value Tmin = 2/u; at the exceptional 
point Q = (r/ — > oo). The maximum, Tmax = ^r/w, is obtained for Q = uj. Since Q = u) implies 
rj = 0, this case corresponds to the Hermitian Hamiltonian. Next, we let the energy constraint 
be -E+ — -E- = 0, = const. Then, as follows from equation (4.12), Tp < tt/O,, and, just as above, 
the passage time achieves the maximum Tmax = 7r/J7 at the point lo = Q corresponding to the 
Hermitian Hamiltonian. For uj ^ we obtain Tp ~ 2/uj, and Tp vanishes in the limit oj oo. 

Similar consideration of the distance between the initial and final states yields 

Lp = LOT = arctan 



and, under the constraint uj = const, we have 2 < L < vr. The upper bound -Lmax = tT; 
being identical to the geodesic distance between the initial and final states defined either on the 
Bloch sphere S'^ or on the one-sheeted hyperboloid B[2 [24], is thus achieved for the Hermitian 
Hamiltonian (O = w). The lower bound, Lmin = 2, is obtained at the exceptional point ($7 = 0). 
This agrees with our general conclusions regarding the behavior of the system in the vicinity 
of the exceptional point (see Section 3). Next, imposing the constraint 0, = const, we obtain 
the same result 2 < Lp < vr. The minimum of Lp corresponds to the limit uj ^ 0,, and, just as 
above, the upper bound is achieved for the Hermitian Hamiltonian (uj = Q). 

We now turn our discussion to the recent controversy around the possibility of achieving 
faster evolution in a quantum-mechanical system governed by a non-Hermitian "PT-symmetric 
Hamiltonian as compared to the equivalent Hermitian system [2, 6, 7, 9]. It should be noted that 
the critique of the results obtained in [2] is essentially based on the following theorem [7] : The 
lower bound on the travel time {upper bound on the speed) of unitary evolutions is a universal 
quantity, independent of whether the evolution is generated by a Hermitian or a non-Hermitian 
Hamiltonian. Analyzing the proof of this theorem, one can see that it is founded on the following 
assumption: the minimal travel time is realized by quantum evolution along the geodesic path 
in the Hilbert space joining initial and final states. 
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This is true in the case of the Hermitian Hamiltonian; however, for a non-Hermitian Hamil- 
tonian, the situation is quite different [24]. Let rij and Uf denote antipodal states on the Bloch 
sphere 5^, and rrij and mj denote corresponding (antipodal) states on the one-sheeted hyper- 
boloid H^. Then, the geodesic distance Lg between rij and nj calculated over the Bloch sphere S'^ 
is identical to the geodesic distance between rrij and mj computed over the one-sheeted hyper- 
boloid (for detailed calculations, see [24]). Moreover, under the same set of constraints, the 
amount of time Tg required to evolve rij into on the Bloch sphere is equal to that required 
to evolve rrij into on by geodesic evolution. This is in accordance with the conclusions 
made by Mostafazadeh in [7]. However, as we have shown above, in the case of the Hermitian 
Hamiltonian, Tg is the lower bound on the time evolution, and, for the non-Hermitian Hamilto- 
nian, it yields only the upper bound on the time evolution. Thus, in non-Hermitian quantum 
mechanics, the evolution of a system is indeed faster than in Hermitian quantum mechanics, 
subject to the same energy constraint. 

5 Conclusion 

In this paper, we considered the non-Hermitian quantum brachistochrone problem for the generic 
non-Hermitian Hamiltonian and focused our attention on the geometric aspects of the problem. 
We demonstrated that for the generic non-Hermitian Hamiltonian, the quantum brachistochrone 
problem can be effectively formulated on the complex Bloch sphere 5"^, enabling the latter to 
be considered as the quantum phase space of the related two-level system. In particular, for the 
effective non-Hermitian Hamiltonian with a real energy spectrum, the corresponding quantum 
phase space is represented by the one-sheeted hyperboloid BP. 

As noted in [4, 5], in the Hermitian quantum brachistochrone problem, the lower bound on 
the time evolution interval r, being proportional to the geodesic distance between the initial 
and final states on the Bloch sphere, is determined using the Fubini-Studi metric on CP^ = S'^. 
We have shown that the geodesic distance Lg between antipodal states rrj and nj calculated 
over the Bloch sphere S*^ is identical to the geodesic distance between corresponding antipodal 
states rrij and uij calculated over the one-sheeted hyperboloid H^. Moreover, the amount of 
time Tg required to evolve rr^ into rrj on the Bloch sphere is equal to that required to evolve rrij 
into nif on by the geodesic evolution. This is in accordance with the conclusions made in [7]. 
However, in the case of the Hermitian Hamiltonian, Tg is the lower bound on the time evolution, 
and, for the non-Hermitian Hamiltonian, it yields only the upper bound on the time evolution. 
Furthermore, while for a quantum-mechanical system governed by the Hermitian Hamiltonian 
the evolution, speed v is bounded by f < Vg, where Vg is the speed of the geodesic evolution, for 
a non-Hermitian quantum system with the same energy constraint, we have v > Vg. This proves 
that in non-Hermitian quantum mechanics, evolution can be faster than in Hermitian quantum 
mechanics [2]. 
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